 [p0,err,k,y] = Newton('f','df',-1,1.0e-6,1.0e-6,1000)